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The variational cluster approximation proposed by Potthoff is applied to the calculation of the 
single-particle spectral function of the transition metal oxides MnO, CoO and NiO. Trial self-energies 
and the numerical value of the Luttinger-Ward functional are obtained by exact diagonalization of 
a TM06-cluster. The single-particle parameters of this cluster serve as variational parameters to 
construct a stationary point of the grand potential of the lattice system. The stationary point 
is found by a crossover procedure which allows to go continuously from an array of disconnected 
clusters to the lattice system. The self-energy is found to contain irrelevant degrees of freedom which 
have marginal impact on the grand potential and which need to be excluded to obtain meaningful 
results. The obtained spectral functions are in good agreement with experimental data. 
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I. INTRODUCTION 

The theoretical description of compounds containing 
partially filled 3d, 4/ or 5/ shells is a much-studied prob- 
lem in solid-state theory. Due to the small spatial extent 
of these shells the Coulomb repulsion between the elec- 
trons in the conduction bands formed from these shells 
becomes unusually strong and approximations which rely 
on a mapping of the physical electron system onto one of 
fictious free particles in a suitably constructed effective 
potential - as is the case in density functional theory [l| 
in the local density approximation (LDA) - cannot even 
qualitatively describe the resulting state. A frequently 
cited example are the transition metal (TM) oxides NiO, 
CoO and MnO. Band structure calculations for the para- 
magnetic phase predict these materials to be metallic 
while experimentally they remain insulators well above 
their respective Neel temperature. For NiO and MnO an 
insulating ground state can be obtained in the framework 
of band theory by introducing antiferromagnetic order - 
for CoO on the other hand even the antiferromagnetic 
ground state is metallic 0- From a comparison of X-ray 
photoemission spectroscopy (XPS) and bremsstrahlung 
isochromat spectroscopy (BIS) it was found[| that the 
band gap predicted by LDA for antiferromagnetic is too 
small by a factor of « 10. Moreover electron spectroscopy 
shows that the electronic structure remains essentially 
unchanged at the Neel temperature for both NiO[i[ and 
CoO@. In both compounds there is practically no dif- 
ference between the electronic spectra in the antiferro- 
magnetic and paramagnetic phase. The same holds true 
for the related compound NiS where LDA band structure 
calculations on the contrary predict that the transition 
to the magnetically ordered phase is accompanied by a 
significant change of the electronic structure p|. 
Failure to predict the insulating ground state and the 
magnitude of the insulating gap is not the only shortcom- 
ing of LDA. In valence band photoemission spectroscopy 
(PES) all three oxides NiO, CoO and MnO show a 'satel- 
lite' at an energy of « 6 — 8eV below the valence band 
top[E 0, The Fano-like intensity variation with pho- 



ton energy at the TM 3p — ► 3d threshold identifies this 
feature as being due to d n — > d"" 1 transitions [9]. This 
part of the electronic structure is not at all reproduced 
by band struture calculations which on the contray would 
predict the d n — » transitions near the valence band 
top. And finally experimental band structures measured 
by angle resolved photoelectron spectroscopy (ARPES) 
show that for all compounds, NiO, CoO and MnO, the 
TM-derived bands near the valence band top are almost 
dispersionless[E H, [nj. This is also in contradition to 
LDA calculations which predicts band widths of around 
2eV for the TM3d-derived bands. 

Starting with the work of Hubbard[ll| a variety of the- 
oretical methods have been invented to deal with this 
problem[ll [H, Q EE M, E3, M M EE EH El . Major 
progress towards a quantitative description of 3d TM ox- 
ides has been made by the cluster method initiated by Fu- 
jimori and MinamiJiEEIEEEEEl This takes the op ~ 
posite point of view as compared to band theory, namely 
to abandon translational invariance and instead treat ex- 
actly - by means of atomic multiplet theory [IE Eil - the 
Coulomb interaction in the 3d-shell of a TM-ion in an 
octahedral 'cage' of nearest-neighbor oxygen atoms. The 
angle-integrated valence band photoemission spectra cal- 
culated by this method are in excellent agreement with 
experiment [H EI EE EE E3- This is clear evidence 
that the atomic multiplets of the partly filled 3d-shell - 
suitably modified by the crystalline electric field (CEF) 
- persist in the solid and play an important part in the 
physics. On the other hand due to its 'impurity' char- 
acter the cluster method can only give k-independent 
quantities. 

Recently, however, ideas have been put forward to 
broaden the correlated ionization and affinity states of 
finite clusters into bands [3fl l3ll. |32| . A particularly ele- 
gant way to do so - the variational cluster approximation 
(VCA) - has been proposed by Potthoff [33j| . Building 
on field-theoretical work of Luttinger and Ward 34( who 
showed that the grand canonical potential f2 of an inter- 
acting Fermion system is stationary with respect to vari- 
ations of the electronic self-energy £(w), Potthoff pro- 
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posed to generate trial self-energies numerically by ex- 
act diagonalization of finite clusters and use them in a 
variational scheme for f2. This amounts to finding the 
best approximation to the true self-energy of the lattice 
amongst the subset of 'cluster representable' ones, i.e. 
exact self-energies of finite clusters. 
So far the VCA has been applied mainly to simplified sys- 
tems such as the single-band Hubbard- model [33J, 135J but 
the success of the cluster method for TM-oxides clearly 
suggests to apply the VCA also to a realistic model for 
TM-oxides thereby using the octahedral clusters intro- 
duced by Fujimori and Minami to generate self-energies. 
Here we outline such a calculation for NiO, CoO and 
MnO, which all have the rocksalt structure. For simplic- 
ity we neglect any of the lattice distortions observed in 
the actual compounds as well as the antiferromagnetic or- 
der at low temperature and study the ideal rocksalt struc- 
ture in the paramagnetic phase. As already mentioned 
the single-particle spectra of these compounds do not 
change appreciably during the Neel transition so this is 
probably a reasonable assumption. A preliminary st udy 
for NiO using the VCA has been published elsewhere |36j. 
Using clusters containing just a single TM-ion implies 
that the self-energy is site-diagonal, i.e. fc-independent. 
The corresponding approximation thus is similar to the 
dynamical mean-field (DMFT) calculations which have 
recently been applied to a variety of compounds [37| . The 
relationship between DMFT and the VCA has been dis- 
cussed in detail by Potthoff[33| and a detailed comparison 
with recent DMFT calculations for NiO will be presented 
below. 



II. HAMILTONIAN 

We start by defining the Hamiltonian which describes 
the correlated TMO lattice. We use a linear combination 
of atomic orbitals (LCAO) parameterization of the non- 
interacting Hamiltonian with hybridization integrals ob- 
tained from a fit to a LDA band structure thereby closely 
following the procedure outlined by Matthciss 38]. The 
parameters so obtained are listed in Table |TJ To give an 
impression about what accuracy can be expected from 
such a fit Figure Q] shows the actual LDA band structure 
and the LC AO-fit for CoO. Following Mattheiss[n| an 
02s orbital was included into the basis set in addition to 
the 02p and TM3d orbitals. This turned out to be cru- 
cial for a correct fit of the dispersion of some TAf3d-like 
bands along T — X. The energy of the 02s orbital has no 
particular impact on the dispersion of the bands near the 
Fermi level and was kept at lAeV below the 02p energy. 
The fit can be improved substantially by including also 
nonvanishing overlap integrals between 02p and TM3d 
orbitals - since, however, the VCA needs well defined 
TMZd orbitals to which the self-energy can be added 
these overlap integrals were omitted. The hybridization 
element (ddS) turned out to be meaningful only in combi- 
nation with these overlap integrals - a fit without overlap 





NiO 


CoO 


MnO 


(ppa) 


0.695 


0.627 


0.542 


(ppn) 


-0.118 


-0.111 


-0.108 


(sda) 


-1.210 


-1.210 


-1.319 


(pda) 


-1.289 


-1.276 


-1.275 


(pdn) 


0.614 


0.596 


0.587 


(dda) 


-0.255 


-0.274 


-0.331 


(ddiv) 


0.060 


0.067 


0.097 


e(02s) 


-14.000 


-14.000 


-14.000 


e(02p) 


0.000 


0.000 


0.000 


e(TM3d) 


2.822 


3.400 


3.899 


WDq 


0.138 


0.142 


0.069 



TABLE I: Hybridization integrals and site-energies e (in eV) 
obtained by a LCAO fit to paramagnetic LDA band struc- 
tures. The site-energies have been shifted so as to have 
e(02p) = 0. 



produced a positive value of (ddS) - and hence was set to 
be zero. By and large the variation of the hybridization 
integrals along the series NiO^MnO is consistent with 
the increasing lattice constant on one hand and the in- 
creasing d-shell radius on the other. The parameters are 
similar to previous estimates in the literature [23l l24j . 
The Coulomb interaction within the d-shell can be writ- 
ten as 

#i= E K% 4 d iA 2 d K3 d Ki . (i) 

Here we have suppressed the site label i and n — (a, a) 
where a G {d xy ,d xz . . .ds z 2_ r 2} denotes the CEF-level. 
The matrix elements can be expressed [H, in 

terms of the 3 Racah-parameters, A, B and C. Due 
to the 'breathing' of the 3d radial wave function these 
parameters should be taken to depend on the d-shell 
occupation[23]. This, however, would create an 'implicit' 
interaction containing terms higher than quartic in the 
Fermion operators and also an interaction between d and 
p-electrons. This would defeat our formalism and we 
therefore do not take the dependence on <i-shell occu- 
pation into account. Whereas B and C can be estimated 
from atomic Hartree-Fock wave functions the parameter 
A is reduced substantially from its atomic value by solid 
state screening. This parameter can in principle be ob- 
tained from density functional calculation s I39L 1401. I n the 
cluster calculations for TM oxides [H [U [gl]2i||27[ A is 
usually treated as an adjustable parameter and here we 
do the same. Adjusting A is equivalent to adjusting the 
'Hubbard [/' = E(d n+1 )+E(d n - 1 )-2E{d n ) where E(d n ) 
denotes the energy of a (i-shell with n electrons. There is 
some ambiguity as to what exactly is to be understood by 
'the energy of c?™' - here we follow Refs. [24j, |25|, |26[ and 
take E(d n ) to be the Coulomb energy of the Hund's rule 
ground state of the free ion i.e. calculated without CEF 
splitting. E(d n ) - and hence U - then can be expressed in 
terms of the Racah parameters (29j .Another way to define 
a Hubbard U would be to note that the average Coulomb 
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FIG. 1: LDA band structure (top) and LCAO-fit (bottom) 
for CoO. 

energy of the d n multiplets is[H, [29j 
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which would suggest to define U av — A — -g-B + |C. 
A second parameter which of importance for charge 
transfer systems[4lj which is usually adjusted to exper- 
iment is the d-level energy or equivalently the charge 
transfer energy A = E(d n+1 L) — E(d n ). Expressing 
the E(d n ) in terms of Racah-parameters A can be ex- 
pressed in terms of these and the difference of site ener- 
gies e(TM3d) — e(02p). The values used in the present 
calculation are given in Table HJ as well. It should be 
noted that while the LCAO-fit actually gives an en- 
ergy for these site energies - see Table U - these values 
of e(TM3d) already incorporate the Coulomb interac- 
tion between d-electrons. Since in our formalism the 
Coulomb interaction is described by the Hamiltonian ((T|) 
one would have to subtract it off to avoid double count- 
ing. Namely if one considers the Hubbard U as known 
one may estimate the 'bare' value of e(TM3d) from the 
LDA site-energy e(TM3d)LDA and the d-elevel occupan- 



cies rid as[20j 

i d = e{TM3d) LDA -9U av n d . (3) 

The estimates obtained in this way are also given in Ta- 
ble im For NiO and CoO these corrected LDA values 
are close to the adjusted parameters used in the actual 
calculation. The situation is different for MnO but for 
this compond the U obtained from ground state energies 
also differs strongly from the multiplet average U av . The 
reason is the strong exchange stabilization in the high- 
spin ground state of d 5 and one may not hope to obtain 
agreement between the two estimates for the site energy 
either. The reason is simply that a 'Hubbard V is not 
uniquely defined in the presence of strong multiplet split- 
ting. Finally, any Coulomb interaction between electrons 





NiO 


CoO 


MnO 


A 


8.25 


7.2 


6.1 


B 


0.13 


0.14 


0.12 


C 


0.60 


0.54 


0.41 


e{TM3d) 


-62.0 


-45.5 


-23.3 


U 


8.38 


7.34 


10.65 


A 


7.42 


11.48 


10.07 


Uav 


8.51 


7.40 


6.23 


nd(LDA) 


0.85 


0.76 


0.56 




-60.53 


-45.81 


-26.91 



TABLE II: Racah parameters, d-level energy, Hubbard U and 
charge transfer energy A. Also given are the 'average Hub- 
bard U' according to [2] the electron number per d-orbital as 
obtained from the LDA calculation and the estimate for the 
d-lebel energy according to [3] All energies are in eV. 

which arc not in the same TM3rf-shell is neglected. 

III. VARIATIONAL CLUSTER 
APPROXIMATION 

Having specified the stronly correlated problem under 
discussion we outline the variational cluster approxima- 
tion. This is based on an expression for the grand poten- 
tial of an interacting many-Fermion system derived by 
Luttinger and Ward 34] . In a multi-band system where 
the Green's function G(k,a->), the noninteracting kinetic 
energy t(k) and the self-energy S(k, ui) for given energy 
lo and momentum k are matrices of dimension 2n x 2n, 
with n the number of orbitals in the unit cell, it reads [13] 

n = -\ e^° + ln det (-G^pc, w„) + (4) 

where lo v = (2v + l)ir//3 with (3 the inverse temperature 
are the Fermionic Matsubara frequencies, 

G- 1 (k,w)=w + A t-t(k)-S(k,w). (5) 

with \i the chemical potential and the functional F[S\ 
is the Legendre transform of the Luttinger- Ward func- 
tional $[G]. The latter is defined [341] as the sum of all 
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closed linked skeleton diagrams with the non-interacting 
Green's functions replaced by the full Green's func- 
tions. A nonperturbative derivation of a functional with 
the same properties as $ has recently been given by 
Potthoff |43j . 3>[G] is the generating functional of the 
self-energy X i.e. 



E « 



(k, w v ) 



9$ 



dGji{k,oj v y 



(6) 



and F[S\ is obtained by Legendre-transform to eliminate 
G in favour of X: 

F[n] = $[G] -\Y. trace (( G ( k ' ^) s ( k > ^)) ■ 



By virtue of being a Legendre transform it obeys 

dF 



G tj (k,u} u ) 



dT,ji(k,LJ u ) 



and using the identity 







dA, 



In det A = (A 



(7) 



(8) 



which holds for any matrix A with det A ^ as well as 
the Dyson equation ([5]) we find that f2 is stationary with 
respect to variations of X: 



0. 



(9) 



The crucial obstacle in exploiting this stationarity prop- 
erty in a variational scheme for the self-energy X is the 
evaluation of the functional F[S\ for a given 'trial X'. 
Potthoff's solution [331] makes use of the fact that just 
like 3>[G], F[E] has no explicit dependence on the single- 
particle terms of H and therefore is the same functional 
of X for any two systems with the same interaction part 
of the Hamiltonian. This is easily seen from the diagram- 
matic representation of <&[G] because the expression as- 
sociated with a given Feynman diagram involves only the 
interaction matrix elements and the Green's function it- 
self. 

In the VCA this independence of F[X] on the single- 
particle terms of the Hamiltonian is used to construct 
trial self-energies by exact diagonalization of finite clus- 
ters and thereby obtain the exact numerical value of 
F[X] . In a first step one chooses a so-called reference sys- 
tem, which has the same interaction part as the lattice 
problem under study but consists of disconnected finite 
clusters. If the interaction terms are short ranged - which 
is the reason for keeping only the Coulomb interaction 
between electrons in the same cf-shell - this can can al- 
ways be achieved by suitable choice of the single-particle 
terms. The disconneted finite clusters of the reference 
system then are solved by exact diagonalization, which 
gives the eigenenergies e„ and wave functions |$„) for all 
particle numbers in the cluster. Of course this sets some 



limit on the size of the clusters. Next, the Green's func- 
tion G(cj) and grand potential Q of the reference system 
are calculated numerically and equation (g]) is reverted 
to express the exact numerical value of F[X] in terms 
of these. This simply means that the summation of in- 
finitely many Feynman diagrams and Legendre transform 
is done implicitely in the course of the exact diagonaliza- 
tion of the reference system. Then, the self-energy X(oj) 
of the reference system - which is readily extracted from 
the Dyson equation for G(w) - can be used as a trial 
self-energy for the lattice system. Thereby the numerical 
value of F[X] calculated in the cluster is simply inserted 
into the Luttinger-Ward formula([4]) for the grand poten- 
tial of the physical (i.e. lattice) system. The variation of 
X(w) is performed by varying the single-electron param- 
eters - such as hybridization integrals or site-energies - of 
the reference system. 

In applying this procedure one frequently has to evaluate 
expressions of the type (the momentum k is suppressed 
for brevity) 



S = 



e^ 0+ In det (-GT 



L M)- (io) 



To evaluate this we closely follow Luttinger and Ward [34 
and first convert the sum into a contour integral 



f(u)g(u)du (11) 



where /(w) is the Fermi function and Co is the standard 
contour encircling the singularities of the Fermi function 
in counter-clock-wise direction. Next we deform the con- 
tour so that it encircles the singularities of the logarithm, 
which arc all located on the real axis (see the Appendix) . 
Following Luttinger and Ward we then use 



(12) 



and integrate by parts. Using (jHJ we thus obtain 
1 



2nf3i 



duo log(l + e^) 

trace [(1 - ^^-)G(u)] (13) 
duj 



where C is a contour that encircles the singularities of 
trace in clockwise fashion. This can now be evaluated by 
numerical contour integration. To derive the expression 
given by Potthoff we note the alternative expression 



S 



-1 
2n(3i 



duj log(l 



c 



)E 



1 dhjuj) 
Xi(uj) duj 



where Xi(oj) are the eigenvalues of G{uS), There are two 
types of singularities of this expression: 
a) zeroes of an eigenvalue, i.e. 
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FIG. 2: Integration contour for the numerical evaluation of 
S. The dots on the imaginary axis denote the Matsubara 
frequencies, the dotted lines are branch cuts of log(l + e _/3 "). 
Crosses denote poles of the Greens function. 



b) singularities of an eigenvalue, i.e. 

A(w) w _^ > _J_<9A(w) = 1_ 

uj — rfu X(lu) duj lu - rift 

Whence 

s = -\ (e ^ + - E ^ + e ~^)) ' 

i.e. the expression derived by Potthoff [33j . 

For the numerical evaluation of S a slightly different 
procedure is more convenient. For Matsubara frequen- 
cies ui v with \v\ < v max the respective terms in the sum 
(|10[) are evaluated directly by computing the eigenvalues 
of G~ x (uj). For \v\ > v m ax we switch to a contour in- 
tegral using and deform the integration contour as 
indicated in Figure [U Along C\ and C4 the integral is 
evaluated numerically again by calculating the eigenval- 
ues of G _1 (w). Along the positive real axis the integrand 
thereby is cut off by the Fermi function. Along C 2 and 
C3 we integrate by parts using again (TT^J) and deform 
the two pieces into the short piece C5 - which is possible 
because the contour encloses no more singularities of the 
integrand. It is not possible to integrate by parts along 
C\ and C4 because these contours cross the lines - indi- 
cated by dashed lines in Figure [2]- where log(l + e~^ w ) 
has branch cuts. The advantage of this procedure is that 
the resulting integration contour C\ — C5 — C4 can be kept 
far from the singularities of the Greens function and self- 
energy on the real axis so that the integrand will always 
be a smooth function and a numerical integration with 
relatively few mesh points (of order 10 3 ) gives accurate 
results. 




FIG. 3: Physical system (top) and reference system (bottom). 
Arrows indicate nonvanishing hybridization. While the phys- 
ical system is a true lattice, the reference system is an array 
of disconnected clusters. The Coulomb interaction between 
electrons in the TM d-shell is the same for both systems, 
however. 



IV. REFERENCE SYSTEM 

Given the excellent results obtained by the cluster 
method for angle-integrated spectraJH, [2J, [2jJ [2(| [2?| 
it seems natural to use clusters which are equivalent to 
TM06 octahedra as reference system. More precisely we 
choose a reference system where each TM3d orbital d a 
is coupled to one 'ligand' orbital L a with these ligands 
in turn decoupled from each other and the interaction 
within the (i-shell given by fTJ). The reference system 
thus is equivalent to an array of non-overlapping iden- 
tical TMOg clusters where each ligand L a corresponds 
to the unique linear combination of 02p orbitals on the 
six nearest O neighbors of a given TM atom which hy- 
bridizes with the TM3d a orbital - see Figure [31 We write 
the single-particle Hamiltonian for a TML5 cluster as 

H single = V(a)(4 tta L ata +H.c.) 

+ J2(E(a) 4,<A, CT + e(a) L^ a (J L a a ) . 

(14) 

The variational parameters thus are the hopping matrix 
elements (V a ), the ligand energies e(a) and the ei-level 
energies E(a) with a £ {e g , t 2g } - in total we thus have 6 
parameters. It is convenient to rewrite the site-energies 
in the reference system as follows 

E(e g ) = e - e x + 3e 2 /5 + e d , (15) 

E{t 2g ) = e ~e 1 -2e 2 /5 + e d , (16) 

e{e g ) = e + ei + 3e 3 /5, (17) 

e(t 2g ) = e + £i-2e 3 /5. (18) 
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FIG. 4: 'Hybrid system' used to find the stationary point of 
ft. 



As shown by Aichhorn et al. [441] optimization of the 'cen- 
ter of gravity' £o ensures that the electron number ob- 
tained by differentiating Q with respect to fj, is equal to 
the result obtained by integrating the spectral function. 
The search for the stationary point of a function of 6 
variables Xi is a difficult task - even more so because the 
stationary point is not a global minimum or maximum 
and is in fact not even a local extremum but a saddle 
point (see below). This problem has motivated the search 
for functionals other than f2[£] which take an extremum 
value at the physical self-energy [45(. We can solve this 
problem in the following way, however: if we have a set 
of parameters which is sufficiently close to the station- 
ary point we can evaluate the derivatives dfl/dXi and 
d 2 il/d\id\j numerically and use the Newton method to 
find the point where dfl/dXi = 0. Next, instead of the 
true lattice system we choose our 'physical system' to be 
a 'hybrid system' which contains both, the (92p-lattice 
and the Ligands for each TM-ion see Figure 0] We take 
the TM - O and TM - TM hybridization to be multi- 
plied by oti, the TM — L hybridization by oli- For at± = 
and oti = 1 we therefore have the reference system itself 
plus a decoupled 02p-lattice. For this system the exact 
stationary point is known - namely the parameters of the 
reference system itself. On the other hand, for a± — 1 
and oi2 = we have the physical lattice system plus the 
decoupled and hence irrelevant ligands and a solution to 
this system is a solution to the lattice system itself. In 
this way we can go continuously from an exactly solvable 
system to the physical lattice system. In practice we vary 
the parameters a±, ai in steps of 0.1 and start the New- 
ton method using the stationary values of the preceding 
step as initial values. When combined with a simple ex- 
trapolation procedure this yields the stationary point of 
the lattice in « 10 steps with < 2 Newton-iterations in 
each step. Obviously such a crossover procedure can be 
formulated for other applications of the VCA as well. 

We next discuss a second important technical point. 
Having found the stationary point we can calculate the 
matrix of second derivatives d 2 Sl/dXidXj and diagonal- 
ize it. Figure [5] shows a scan of through the station- 
ary point of CoO along the principal axes so obtained. 
There are two important things to recognize. First, the 
stationary point is a saddle point but the above crossover 
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FIG. 5: Scans through the stationary point of CoO. The 
normalized eigenvectors Vj of d^ft/dXidXj were muliplied by 
Xi. 



procedure had no problems to find it. Second there are 
certain directions in parameter space where f2 shows only 
an extremely weak variation. This turned out to be true 
in all other cases studied as well. This weak variation 
may either stem from a near-invariance of the self-energy 
under changes of the cluster-parameters or the change 
of the self-energy is appreciable but irrelevant in that 
it does not change the lattice O. The presence of such 
'nearly invariant lines' in parameter space clearly is un- 
desirable in that it may induce numerical instabilities. 
It may happen that small changes of f2 due to e.g. a 
slightly wrong LCAO band structure or even numerical 
inaccuracies may drive the stationary point along these 
lines in parameter space to compensate for them. To 
simplify matters the number of parameters therefore was 
reduced. Inspection of the eigenvectors associated with 
the 'nearly invariant lines' showed, that these were pre- 
dominantly combinations of the parameters e± and £3 in 
(fT8|) . These parameters have almost no influence on 
and hence were not subject to variation. The parameter 
£1 was set equal to zero. In the cluster calculation the 
value for £3 would be 2(ppa) — 2 (ppn) 24} - for simplicity 
the value e 3 = lAeV was used for all three compounds. 
Including ei and £3 into the set of parameters to be opti- 
mized actually turned out to give unsatisfactory results 
for the single particle spectrum. e\ tended to take on 
large positive values whereas £3 usually took large nega- 
tive values. As a net effect this produced spurious pho- 
toemission peaks with very small spectral weight which 
were split off by one or two eV from the remainder of the 
photoemission spectrum resulting in too small gaps and 
poor agreement with experiment. Clearly this is a feature 
of the variational cluster approximation which needs to 
be clarified. It should be noted that reducing the num- 
ber of parameters which are optimized simply amounts 
to restricting the space of trial self-energies. Since op- 
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FIG. 6: Grand potential Q(T) obtained by the VCA for CoO. 



timization of these parameters hardly changes this is 
similar to restricting the degrees of freedom in a trial 
wave-function to the most relevant ones. It then seems 
that an 'overoptimization' of parameters leads to poor 
results but on the other hand the inclusion of irrelevant 
degrees of freedom into a variational wave function may 
also be detrimental for properties of the wave function 
other than the ground state energy. 
Finally, NiO turned out to be a special case. Since the 
ground state of d s in cubic symmetry is ^fg e g the hopping 
integral Vfeg) has practically no impact on £1 because it 
connects filled orbitals. In fact, derivatives of £1 with re- 
spect to V(t2g) turned out to be of order 10 -10 , i.e. well 
beyond the numerical accuracy of the whole procedure. 
V(t2 9 ) was therefore kept at 2(pdir) which again is the 
value expected in the cluster calculation. In a previous 
VCA-study of NiO[36| a different approach was chosen. 
There the hopping integral V(t2 g ) was set set equal to 
zero. This implies that the t2 9 -like ligands are irrelevant 
alltogether and can be discarded from the reference sys- 
tem, so that also the paramaters E{t2 g ) and efog) play 
no more role. Although slightly different LCAO- and 
Racah-parameters were used in this calculation the re- 
sults obtained in Ref. [36[ are very similar to the ones in 
the present study, in particular the bands in the valence 
band top are essentially identical. 



V. RESULTS 

The parameters at the stationary point of have prac- 
tically no dependence on temperature. Figure [6] shows 
the temperature dependence of f2 for CoO. This can be 
fitted very well by fi(T) = tt - fc B Tlog(12). The sec- 
ond term thereby is the entropy due to the degeneracy 
of the A T\ g ground state of d 7 in cubic symmetry. In a 
system with a wide gap this is the expected behaviour of 
il. This is clearly a trivial result but it should be noted 
that for the discussion of a phase transition to a magnet- 
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FIG. 7: Imaginary part of the Green's function (top) and self- 
energy (bottom) for the reference C0O6 cluster. The imagi- 
nary part of the frequency is —O.leV. Parameter values cor- 
respond to the stationary point at 150K. 

ically or orbitally ordered state the correct description of 
this entropy is important because this competes with the 
energy gain due to ordering. 

Next we consider the resulting self-energy. Figure [7| 
shows the spectral density of the C0O6 cluster and the 
self-energy at the stationary point for T = 150 Kelvin. 
Due to the cubic symmetry of the cluster only the di- 
agonal elements of the self-energy are non- vanishing and 
these are identical between all e g and t 2g orbitals, respec- 
tively. Luttinger has shown [Z(| that the self energy has 
a spectral representation of the form 

£(a,) =?7 +£-^- (19) 

where the real matrix rj is actually the Hartree-Fock po- 
tential and the poles C, v are all on the real axis. The spec- 
tral density of the cluster has a well-defined gap around 
u> = between a charge-transfer peak and the upper Hub- 
bard band. The self-energy for both e g and t 2g electrons 
has a strong central peak (indicating a pole Q v with large 
residuum S„) in this gap. Using the Dyson-equation it 
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is easy to convince oneself that such a strong peak in the 
self-energy indeed 'pushes open' a gap in the pole struc- 
ture of the Green's function. Several other prominent 
peaks create additional gaps in the spectral density and 
thus split off the satellite below —9eV. In addition there 
are many small peaks near the top of the valence band. 
Since the poles of the Green's function are 'sandwiched' 
between the poles of the self-energy we thus expect a 
large number of 3d-derived bands with very small dis- 
persion in this energy range. 

Next, we proceed to a comparison of G(k, to) to 
experiment and begin with NiO. Figure [5] compares k- 
integrated spectral densities at T = 150 Kelvin to an- 
gle integrated valence band photoemission spectra taken 
by Oh et aZ.@ at two different photon energies. It is 
known 47] that with decreasing photon energy the inten- 
sity of 02p derived states increases relative to that of 
TMid derived states - the change of the spectra with 
photon energy thus allows to draw conclusions about 
the character of individual peaks. Moreover, final states 
with d n character are enhanced at a photon energy 



FIG. 9: Top: Dispersion of experimental bands measured by 
ARPESd_C§ in NiO. Bands in the satellite region are given 
with error bars due to their strong broadening. Bottom: k- 
dependent spectral function for momenta along F — X in NiO. 
Lorentzian broadening 0.05 eV, d-like weight is multiplied by 
factor of 2. 



just above the TMip — > TM3d absorption threshold so 
that such energies are particularly suited to identify this 
type of final state. Accordingly, at hv — I50eF the ex- 
perimental spectrum mostly resembles the d-like spectral 
density, whereas at hv = 67eV the states at the valence 
band top are anti-resonantly supressed - whence 02p- 
derived features become more clearly visible - and the 
'satellite' at — lOeV is resonantly enhanced. Figures [S] 
and 1101 compare the fc-resolved spectral function for mo- 
menta along (100) (r -> X) and (110) (r K) to the 
experimental band dispersion by Shen et a?.[ic|]. 
The spectral density has gap of approximately 4eV" 
around the chemical potential. This is consistent with 
experiment [3| but has of course been achieved by the 
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FIG. 10: Same as Figure [9] but for momenta along the (110) 
direction (r — K). 



choice of A and A. At the top of the photoemission 
spectrum, E < 0, there is a high-intensity band complex 
at binding energies between « —3.5 eV and ~ —2 eV, 
which was shown to consist of several sub-peaks by Shen 
et aL[l(|. These authors did not actually resolve the 
dispersion of the individual sub-peaks although the data 
seem to indicate a weak overall 'upward' dispersion as 
one moves r — > X which would be consistent with the- 
ory. Proceeding to more negative binding energy the ex- 
perimental band structure shows a gap of ~ leV and 
then a group of dispersionless bands between —6eV and 
—AeV. This is bounded from below by a weakly disper- 
sive band which resembles one of the 02p derived bands. 
In the angle-integrated spectrum, Figure [5k., the topmost 
of these dispersionless bands produces the 'shoulder' at 
—AeV. The gap between the topmost band complex and 
the group of dispersionless bands in the theoretical spec- 
tra is not as wide as in experiment but there are clearly 



several dispersionless bands in approximately the right 
energy range. The agreement would be very good if the 
peaks at w — 3.5eV in Figures [9] and \W\ were shifted 
by w O.beV to more negative binding energy. In ad- 
dition, the 02p derived band can be seen clearly. As 
can be seen from the band structure in Figure [1] this 
band actually has a saddle point at X - this gives rise to 
a van Hove-singularity in the angle-integrated spectrum 
which matches very well the peak at ~ — 5.5eV in Figure 
[SJd. The sole strongly dispersive feature in the spectrum, 
namely an 02p-derived band at binding energies between 
—6eV — > —9eV is again well reproduced by theory. Fi- 
nally the 'satellite' at binding energies — 8eV — > —12eV 
consists of at least two sub-peaks as can be seen in Figure 
|5b and also in the ARPES data. 

Next we consider CoO. Figure [TT1 compares the angle 
integrated spectra at different photon energies and the k- 
integrated spectral function obatined by the VCA, Figure 
[T2l shows the dispersion along r — ► X and ARPES data 
from Shen et aLQ and Brookes et al.\4^. The XPS spec- 
trum for CoO starts out with a prominent peak at —3eV 
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FIG. 12: Top: Experimental band structure along (100) as 
seen in ARPES by Shen et aZ[f| and by Brookes et o/[4g| in 
CoO. Bottom: fe-dependent spectral function for momenta 
along r — X for CoO. Lorentzian broadening 0.05 eV, d-like 
weight is multiplied by factor of 2. 



followed by three 'humps' at ~5eV, — 8eV and —\2eV. 
The VCA gives peaks of d-like spectral weight at roughly 
these energies although the peak at —5eV is at slightly 
too negative energy. The PES spectrum at 40eV shows 
additional peaks at —6.5eV and —9eV which were in- 
terpreted as 02p-derived by Shen et a2.p|. These peaks 
are also reproduced by the VCA. A little more prob- 
lematic is the angle resolved spectrum. Along (100) the 
VCA predicts a split peak at the the top of the valence 
band at —AeV. This splitting is not seen in experiment 
- on the other hand, the spectra were taken at low pho- 
ton energy where the Co3d states have relatively small 
weight. There is another <i-derived band at —6eV which 
corresponds to the second 'hump' in the angle integrated 
spectrum. This is crossed by and mixes with one of the 



02p derived bands, which start at T at an energy of 
— 5eV. The presence of these two crossing bands may 
explain the 'wiggly' nature of the bands observed exper- 
imentally in this energy range. The presence of more 
than one band and a possible crossing between these is 
clearly seen in the data of Shen et ai. Finally there is 
the strongly dispersive 02p derived band at energies of 
around —8eV. Surprisingly the experimental dispersions 
for this band differ somewhat - this may be due to the 
crossing of this band with the dispersionless Co3d de- 
rived band at —7.5eV which leads to a hybridization gap 
in the 02p derived band. 

Finally we consider MnO. Figure [13] shows the angle 
integrated photoemisison spectra compared to the result 
from the VCA. At high photon energy the experimental 
spectrum matches well the d-derived density of states. 
The single-particle gap and the structure of the valence 
band spectrum are reproduced well. At a photon energy 
of 20eV the intense peak at —5.5eV almost disappears 
and another large peak at — 6.5eV appears, which accord- 
ingly must have 02p character. In the theoretical spectra 
this is reproduced well, the peak at — 6.5eV again is due 
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FIG. 14: Spectral function for momenta along F — X in MnO. 
Lorentzian broadening 0.05 eV, d-like weight is multiplied by 
factor of 2. There are no ARPES results available. 
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to a van-Hove singularity at X. Figure [14] shows the 
k-resolved spectrum along (100). Lad and Hcinrich[8] 
performed ARPES measurements on MnO but did not 
perform any band mapping due to the broad nature of 
peaks so that an experimental dispersion unfortunately 
is not available. 

Lastly we discuss the fine structure of the TM3d derived 
bands. This is shown in Figure 1151 The band struc- 
tures of all three compounds have a similar structure: 
at the top of the band structure there is a group of dis- 
persive bands which shows a rough similarity with the 
upper group of bands in the LDA band structure (see 
Figure [1]), which have mainly TM 3d-character. The to- 
tal width of this band complex is reduced by a factor 
of m 0.5 as compared to LDA. These bands have high 
spectral weight and produce the intense peaks at the top 
of the angle- integrated spectra for NiO and CoO. Sepa- 
rated from this group of dispersive bands there is then a 
region with many almost dispersionless bands with rela- 
tively low spectral weight. This overall structure can be 
understood by considering the spectral representation of 
the self energy (jTHJ) and the equation for the poles of the 
Green's function 



co - 



[i — £fc — Re X(w) = 



FIG. 15: Dispersion of the 'sub-peaks' in the photoemission 
spectra of NiO (top) and CoO (bottom). Also shown are the 
energies of the poles of the self-energy. 



where we have considered the single band case for sim- 
plicity. Since Re S(w) takes any value between oo and 
—oo in between two successive poles Q v and C^+i there is 
one band in between any two successive poles of the self 
energy. This implies that the distance between these two 
successive poles is an upper bound for the width of this 
band, which may be viewed as a kind of correlation nar- 
rowing. Moreover, if a pole £ v has only a small residuum 
the resulting pole of the Green's function will be almost 
'pinned' very close to it - as can be seen repeatedly in 
Figure 1151 The topmost group of relatively strongly dis- 
persive bands then is actually above the topmost pole of 
in the valence band region and the dominant 'gap 
opening peak' in the center of the insulating gap, see Fig- 
ure[7] Since the separation in energy between these peaks 
is large - of the order of the insulating gap - these bands 
still have an appreciable width. The similarity with the 
LDA band structure is due to the fact that the dispersion 
of these bands is largely due to direct d- ri-hopping, which 
remains operative also when the self-energy is included. 
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FIG. 16: fc-integrated d-like spectral weight obtained by 
DMFT (top) (from Ref.Q and by VCA (bottom) compared 
to XPS-data from NiO. 



While the fine structure of the valence band top is not 
really resolved experimentally as yet, at least experiment 
puts a quite low upper limit - < 0.5eV along (100) - on 
the width of the individual bands. The VCA would be 
consistent with that. 

The large number of dispersionless bands at more nega- 
tive binding energy is produced by the large number of 
densely spaced poles of the self-energy. These are in turn 
the consequence of the large number of CEF-split multi- 
plet states in the TM06 cluster. Interestingly, at least in 
the case of NiO where detailed band mapping is available 
from ARPES, the experimental band structure is quite 
consistent with this overall structure, namely a group 
of dispersive bands at the top of the valence band and 
essentially dispersionless bands at more negative bind- 
ing energy. The distance between the dispersive band 
complex and the dispersionless bands is underestimated 
somewhat by VCA. 

To conclude we compare the results of the VCA for NiO 
with recent LDA+DMFT calculations [U HI]. Figure 
[TBI compares the k-integrated spectra, Figure [T7] shows 
the dipersion along (100). The DMFT-results are taken 
from Yin et al. [22| which have obatined essentially iden- 
tical results as Kunes et a/.[21[. While the k-integrated 
spectra look similar at first sight, comparison with the 
band structure shows that there are major differences. 
In the DMFT spectrum the top of the valence band is 
formed by a split off peak A with low spectral weight. 
This corresponds to the two topmost bands labelled A in 
Figure [T7] This form of the density of states is actually 
reminiscent of the results of the three-body-scattering 



theory of Manghi et aL[l5(. In the VCA spectrum the 
top of the valence band is formed by an intense peak A' , 
which corresponds to the topmost band complex A' in 
Figure [T71 The intense peak B in the DMFT spectrum 
on the other hand originates from the band B in Figure 
n~7l The DMFT bands moreover show a rather obvious 
correspondence with the LDA band structure, resulting 
in bands with quite strong dispersion. As already noted 
the VCA differs strongly from LDA and shows a larger 
number of bands, with several of them being practically 
dispersionless, i.e. corresponding to localized electrons. 
Comparing with experiment, the raw data of Shen et 
al. [h} show no indication for the split off bands A with 
low spectral weight as predicted by DMFT. With the 
exception of the 02p derived bands ARPES moreover 
shows no indication of the wide bands predicted by 
DMFT - rather there is a number of dispersionless bands 
as expected on the basis of the VCA. One may say that 
there are major differences between DMFT and VCA so 
that further experiments might resolve this discrepancy. 



VI. CONCLUSION 

To summarize, the variational cluster approximation 
due to Potthoff allows to combine the powerful clus- 
ter configuration interaction method for transition metal 
compounds due to Fujimori and Minami with the field- 
theoretical work of Luttinger and Ward to implement a 
variational scheme for the electronic self-energy and con- 
struct an efficient band structure method for strongly 
correlated electron compounds. As demonstrated above, 
a realistic band structure and the full atomic multiplet in- 
teraction can be incorporated into the Hamiltonian with- 
out problems, the system can be studied at arbitrar- 
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ily low temperatures and the Green's function be ob- 
tained with arbitrary energy resolution. The key numer- 
ical problem of finding the stationary point of the grand 
potential thereby can be solved efficiently by a simple 
crossover procedure. It has been shown that in the course 
of varying the self-energy there may exist redundant de- 
grees of freedom which leave the grand potential almost 
unchanged. Such redundant degrees of freedom can be 
eliminated by simply reducing the numer of parameters 
in the reference system. 

The results are quite encouraging in that there is good 
agreement between the calculated Green's function and 
electron spectroscopies at least to the extent that ARPES 
data are available. The good agreement also suggests 
that the band structure of NiO, CoO and MnO is 
'Coulomb generated' in that the atomic multiplet struc- 
ture survives with minor modifications. All in all the 
VCA appears to be a promising tool for the study of 
realistic models of correlated electron systems. The pos- 
sibility to treat the multiplet and CEF-splitting of the 
various TM3d configurations more or less exactly should 
make it possible to address magnetic or orbital ordering 
phenomena in transition metal compounds. 
I would like to thank M. Potthoff for many instructive 
discussions. 

VII. APPENDIX 

We show that ln(det G(u>)) is analytical off the real 
axis, where G(u>) can be either the exact Greens func- 
tion of the reference system or the approximate Greens 
function form the VCA. It is sufficient to proove that all 
eigenvalues of G(uj) have a nonvanishing imaginary part 
for oj not on the real axis. This is prooved in turn if we 
show that 

(v|G(w)|v) =^«i Gij(u) Vj 



has a nonvanishing imaginary part for any normalized 
v. For the exact Greens function we have - using the 
Lehman representation - 

M°MI->4l ,-^) ''* + '*» 

where e„ = E v — fj,N v and 

C V ' V = (v'\'Y]viCiW) 

i 

It follows that for u> in the upper (lower) half plane all 
eigenvalues of G(uj) have a negative (positive) imaginary 
part and accordingly all eigenvalues of G _1 (w) have a 
positive (negative) imaginary part. The imaginary part 
could only be zero if all C v / V were zero which is unlikely 
to occur. A similar proof has been given previously by 
Dzyaloshinskii[49j. Luttinger has shown that the self- 
energy S(w) has a spectral representation of the form 

with a real g[S|. Since G(lj) is Hermitean for real uj the 
matrices S„ are Hermitean too moreover positive definite. 
Namely if one of the matrices S„ had a negative eigen- 
value A then G -1 ^ + ie) had the eigenvalue ~A plus 
terms which stay finite as e — > 0, whereas we have shown 
that all eigenvalues for u> in the upper half-plane have 
positive imaginary parts. It follows immediately that 

(v|w - S(w)|v) 

has a positive (negative) imaginary part for lu in the up- 
per (lower) half plane which prooves that all eigenvalues 
of the aproximate Greens functions off the real axis have 
nonvanishing imaginary parts as well. 
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